Effective Hydro-dynamic Boundary Conditions for Microtextured Surfaces 
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Understanding the influence of topographic heterogeneities on liquid flows has become an impor- 
tant issue with the development of microfluidic systems, and more generally for the manipulation 
of liquids at the small scale. Most studies of the boundary flow past such surfaces have concerned 
poorly wetting liquids for which the topography acts to generate superhydrophobic slip. Here we 
focus on topographically patterned but chemically homogeneous surfaces, and measure a drag force 
on a sphere approaching a plane decorated with lyophilic microscopic grooves. A significant de- 
crease in the force compared with predicted even for a superhydrophobic surface is observed. To 
quantify the force we use the effective no-slip boundary condition, which is applied at the imaginary 
smooth homogeneous isotropic surface located at an intermediate position between top and bottom 
of grooves. We relate its location to a surface topology by a simple, but accurate analytical formula. 
Since groves represent the most anisotropic surface, our conclusions are valid for any texture, and 
suggest rules for the rational design of topographically patterned surfaces to generate desired drag. 

PACS numbers: 68.08.-p, 68.35. Ct 



Introduction.- The advent of microfiuidics has moti- 
vated the growing interest in understanding and model- 
ing of flows at small scales or in tiny channels. In recent 
years it has become clear that the no-slip boundary con- 
dition at a solid-liquid interface is valid only for smooth 
hydrophilic surfaces [lH^], and for many other systems 
it does not apply when the size of a system is reduced. 
Thus the hydrophobicity of smooth surfaces could induce 
a partial slippage, v s = bdv/dz, where v s is the velocity 
at the wall, b the slip length, and the axis z is normal to 
the surface |5|. This concept is now well supported by 
nanorheology measurements |l|, [2|, |4[ . 

However, only very few solids are molecularly smooth. 
Most of them are rough, often at a micrometer scale. This 
roughness may be induced by some processes of fabrica- 
tion or coating, but microtextures are also found on the 
surfaces of most plants and animals. In particular, many 
solids arc naturally striated by grooves, which can also be 
prepared for specific microfluidic purposes, such as pas- 
sive chaotic mixing [g, |7| . Most studies of flow past rough 
surfaces have concerned poorly wetting liquids for which 
the topography acts to favor the formation of trapped gas 
bubbles (Cassie state) , and to generate superhydrophobic 
slippage [3, |9( . For rough wettable surfaces the situation 
is unclear, and opposite conclusions have been made: one 
is that roughness generates extremely large slip 11011 , and 
one is that it decreases the degree of slippage [llj . Re- 
cent data (supported by simulations |lj|) suggest that 
the description of flow near rough surfaces has to be cor- 




FIG. 1: Sketch of a sphere approaching a model grooved sur- 
face (left) with the example of a typical experimental signal 
and a scanning electron micrograph (taken under an angle) 
of the surface obtained by a soft lithography (right) 




rected, but rather for a separation, but not slip [13 . 
Another suggestion is to combine these two models 

In this Letter we describe how the boundary condi- 
tions can be modified by the surface texture. We focus 
on the case of special interest where this model surface 
is decorated by rectangular microgrooves, i.e. on the sit- 



uation of the largest possible anisotropy of the flow. We 
analyze the hydrodynamic interaction between a smooth 
sphere and a grooved plane, as sketched in Fig. [TJ and 
the texture parameters are systematically varied at the 
micrometer level, in order to investigate their influence 
on a drag force. Our results do not support some previ- 
ous experimental conclusions on a large slip for similar 
systems. Instead, we unambiguously prove the concept 
of an effective no-slip plane shifted down from the top of 
roughness. To the best of our knowledge, this is the first 
study, where experimentally found values of this shift 
were quantified theoretically and related analytically to 
controlled parameters of topographic patterns. 

Experimental.- We use a specially designed home- 
made setup 16l-ll8j to measure on a microscale the dis- 
placement of a sphere towards the corrugated wall at 
constant gravity force. The steel sphere of density p p = 
7.8 x 10 3 kgm~ 3 and radii ranging from 3.5 mm to 6.35 
mm is embedded in a liquid contained in a cylindrical 
glass vessel with a 50 mm diameter and a 40 mm height. 
As a liquid we have chosen high molecular weight PDMS 
(silicone) oil (47V100000 Rhodorsyl oil, from Rhonc- 
Poulenc), with dynamic viscosity fj, = 97.8 Pas at 25°C, 
which is Newtonian for shear rates up to 100 s _1 . Such 
shear rates are never reached in our experiment. 

Table 1 : Parameters of the textured samples and the shift 
of effective hydrodynamic wall, s. 
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The microstructured surfaces were created by com- 
mon soft lithography, in a three steps process, trans- 
ferring geometric shapes from a mask: first to a sili- 
con wafer coated with a (SU8) photoresist, second to a 
replica molding obtained by soft imprint of a thcrmo- 
reticulable PDMS, and finally to a replica of the PDMS 
mold by soft imprint of a thiolene based resin (NOA 81, 
Norland optical adhesives) on glass microscope slides (to 
be fixed at the bottom of the vessel). This resin was 
chosen for its resistance to compression and to solvent- 
swelling, and for its good adhesion to glass |19J. The final 
structures are checked by scanning electron microscopy 
(see Fig. []}. The textures are characterized by spacing 
(5, height e and period L. The liquid fraction, cf> = S/L, 
can be precisely measured since it is the ratio of the up- 
per surface of the crenellations over the total surface of 
the sample. It varies largely with the patterns (from 0.1 



to 0.9), and e/L varies from 0.168 to 0.45, as displayed 
in Table 1. Contact angles against PDMS for all tex- 
tures were found to be below 30° , so that surfaces can be 
considered as lyophilic. Therefore, we expect PDMS to 
invade the surface texture (Wenzel state). 

We measure the distance, h, which is defined from the 
top of the textures (contact) by using an interferometric 
technique [16l4l8| with the accuracy 0.2 /im. The veloc- 
ity U(h) of the sphere is found by multiplying the ve- 
locity of interference fringes displacement by a factor of 
A/2n, where A = 632.8 nm is the wavelength of the He-Nc 
laser, and n = 1.404 the refraction index of PDMS. After 
opto-electronic conversion and amplification, the signal 
is recorded with a high frequency electronic oscilloscope 
(DPO4032 from Tektronics). A deceleration of the sphere 
(Fig. [1]) is reflected in the increase of the period of the 
signal, until contact occurs, and its position is defined 
from the recorded signal, at the time when the period 
of the signal becomes very large indicating a vanishing 
velocity. Note that the signal-to-noise ratio deteriorates 
at vanishing frequency, because the low frequency limit 
of the oscilloscope is reached. The measured frequency is 
averaged over 7 to 8 periods, except just before the con- 
tact, where no averaging is applied in order to capture 
the rapid velocity variations occurring in that region. 

Results and discussion. - Fig. [2] shows the drag F 
(equal to the gravity force) scaled by the Stokes force 
Fst = 6n/j,aU(h), i.e. U(oo)/U(h). Solid line is a the- 
oretical force (Taylor's equation) predicted for a case of 
smooth wall and no slippage at the surface: 



F T /F st = a/h. 



(1) 



Also included are the experimental data for samples with 
similar e, but different <j> and L, which show deviations 
from the behavior predicted by Eq. [1] Close to the wall, 
for a/h > 50, the drag is always significantly less than 
the force near a smooth wall, and this reduction increases 
with <fi. 

To examine these deviations we evaluate a correction 
to the drag force, 



f*(h) = F(h)/F T (h) 



(2) 



Note that in general case for a rough surface /* should 
also depend on the radius of the sphere [l2j. However, 
with our configuration geometry experimental data do 
not vary with a. This is well illustrated in Fig. [3ja), 
where the experimental values of /* obtained with sam- 
ple #7 at different a and plotted as a function of h/L 
collapse into a single curve, which tends to unity at 
large distances and decreases significantly when h be- 
comes of the order of L and smaller. Since at short 
separations we observe /* — > 0, one can conclude that 
slippage (which would lead to /* — > 1/4 [20| ) obvi- 
ously does not mimic roughness when h is small, by 
overestimating the drag force. The same remark con- 
cerns effective supcrhydrophobic slippage where /* — > 
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FIG. 2: Drag (squares) scaled by the Stokes force for walls 
decorated with grooves of similar heights (e = 42 — 45 /im), 
but different cj> and L. From top to bottom the data sets for 
samples #6, #2, #5, and #7 (see Table 1). Solid line shows 
the theoretical prediction for a smooth lyophilic wall, Eq.JT}, 
defined at the top of grooves. The dashed curves from top 
to bottom are calculations of the force expected for a smooth 
lyophilic wall shifted from the textured wall to a distance 
s = 0.5, 4.2, 9 and 28.5 /im. 



2(4 - 30)/(8 + 9(f) - 90 2 ) HI and is equal to ~ 0.3 for 
this particular sample. 

Therefore, our experimental results are now compared 
with theoretical calculations made for an effective smooth 
plane shifted down from the top of the corrugations, i.e. 
by assuming F(h) = Fr(h + s) where s is the value of 
constant, i.e. independent on h, shift. This implies that 



f*(h) 



F T (h + s) _ h 
F T (h) ~ h + s' 



(3) 



Fig. \%1'a) includes a calculation (dashed curve) in which 
an adjustable parameter, a shift of s = 28.5 /im, is incor- 
porated into the Taylor equation. The fit is very good 
for all h, suggesting the validity of the model. Fig. |3jb) 
shows another scries of experiments made with the fixed 
radius of the sphere, but different parameters of the tex- 
ture. If similar fits are made to a variety of experiments 
it is found that the shift of an equivalent plane required 
to fit each run increases from 0.5 /jm for sample #2 to 
28.5 /xm for sample #7. In Table 1 we present the ex- 
perimental values of s for different samples, and theoret- 
ical curves calculated with Eq.© are included in Figs. [5] 
and [3] The fit is excellent at all separations except as 
very small, h/L < 0.01. Thus our experiment shows 
that an effective (scalar) shift, s, is a unique physical pa- 
rameter that fully quantifies drag reduction at a highly 
anisotropic corrugated surface. This striking result indi- 
cates that in our experiment pressure remains isotropic 
despite an anisotropy of the flow. 

Now we try to relate s to parameters of textured sur- 
faces. As proven in 17J, |2l(, for a large gap, h ^> L, the 
shift of the equivalent no-slip plane from the real surface 
is equal to the average of the eigenvalues of the effective 
slip-length tensor (at the slip plane defined at the top of 
asperities) 
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FIG. 3: Measured correction to drag (symbols) (a) for spheres 
of different radii (a = 3.5, 5.75 and 6.35 mm) interacting with 
sample #7. Dashed line shows the calculated correction, using 
s = 28.5 (j,m in Eq.©, (b) for striped walls having nearly same 
heights e = 42 — 45 fim , but different <f> and period L; from 
left to right: samples #6, 2, 5, 7. Dashed lines: model, Eq.©, 
with (a) s - 28.5 pm, (b) from left to right: s = 0.5, 4.2, 9 
and 28.5 \im. 
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FIG. 4: Experimental values of s (symbols) as a function of 
4> for grooves with similar heights (e = 42 — 45 /im), but 
different cf> and L. Lines show theoretical predictions, Eq. 
ffl, where &jj ff and 6^ are calculated by using the analysis of 
Q (solid lines), Eqs. ©-([ill (dash-dotted lines), and Eq. © 
(dashed lines): (a) samples #6 and #7 (black), samples #2, 
#3, (grey); (b) sample #5 (black) and sample #1 (grey). 



Therefore, the problem of calculating s reduces to com- 
puting the two far-field eigenvalues, 6, 



and b^rg, which 



(4) 



attain the maximal and minimal directional slip lengths, 
respectively. 

In the limit e <C L <C h, the theory [23J predicts that 
the effective no-slip surface for arbitrary smooth periodic 



surfaces is at the average height: 



Vfl' 



cj)e, 



(5) 



so that s/e is controlled mainly by </>. To examine the sig- 
nificance of 4> more closely, the experimental s normalized 
by e arc reproduced in Fig. 2] The measured data show 
much smaller s/e than the theoretical prediction of model 
([5]) shown by a dashed line. A possible explanation for 
this discrepancy is that the height of asperities in our ex- 
periments was not small enough, 0.168 < e/L < 0.45. We 
also compared our data with another calculation (solid 
curves) for hydrophilic grooves with finite e/L based on 
numerical results |22| for eigenvalues of the slip-length- 
tensor. Even at moderate e/L theoretical predictions for 
s 22| are much smaller than measured values. 

An alternative model can be obtained if we use the ana- 
lytic solutions for alternating slip and no-slip stripes [24J : 



(0) 




(7) 



where we naturally assumed that the local partial slip 
is equal to the height of grooves. Fig. |4] shows that Eq. 
dU together with Eqs. ((6]) and ([7]) (dashed curves) give 
almost quantitative agreement with experimental data. 
(We also include theoretical values of s to Table 1 to allow 
a direct comparison with experimental results.) There- 
fore, by using the equivalence of a flow past rough and 
heterogeneous surfaces at large scale, we were able to 
quantify a drag reduction at the smaller scale, of the 
order of the size of roughness elements. Note however 
that our results do not apply to a very thin gap situation 
/i«I, where s scales with the channel width [25j, which 
is again consistent with our experiment. 

Concluding remarks.- We have studied a drag force 
on a sphere approaching a corrugated plane. Our ex- 
periment shows quantitatively that in this situation the 
effective no-slip boundary condition, which is applied at 
the imaginary smooth homogeneous isotropic surface lo- 
cated at an intermediate position between top and bot- 
tom of grooves, fully mimics the actual one along the 
true corrugated interface, except as for a very thin gap. 
The location of this effective isotropic plane depends on 
the parameters of the texture, and can be found by us- 
ing simple formulae for effective slip lengths in the limit 
of a thick channel. Since for grooves anisotropy is maxi- 
mized, the same conclusion would be valid for other types 
of anisotropic (e.g., sinusoidal, trapezoidal, and more) 
and/or isotropic (e.g. pillars, etc) textures, but of course, 



Eqs. ([6]) and (JT)) should be replaced by analytical or nu- 
merical solutions for a corresponding texture, as will be 
described in subsequent papers. 

We have also demonstrated that topographically pat- 
terned (Wenzel) surfaces could reduce a drag force more 
efficiently compared to expected even for slipping super- 
hydrophobic (Cassie) textures with trapped gas. There- 
fore, our results suggested rules and a general strategy for 
the rational design of topographically patterned surfaces 
to generate desired low drag. 
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